Theory of (001) surface and bulk states in Yi yCayBa2Cu307_5 
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A self-consistent model is developed for the surface and bulk states of thin 
Yi_j;Caj;Ba2Cu307_5(YCBCO) films. The dispersions of the chain and plane layers are 
modelled by tight-binding bands, and the electronic structure is then calculated for a finite- 
thickness film. The dopant atoms are treated within a virtual crystal approximation. Because 
YCBCO is a polar material, self-consistent treatment of the long range Coulomb interaction leads 
to a transfer of charge between the film surfaces, and to the formation of surface states. The tight 
binding band parameters are constrained by the requirement that the calculated band structure 
of surface states at Cu02-terminated surfaces be in agreement with photoemission experiments. 
The spectral function and density of states are calculated and compared with experiments. Unlike 
the case of Bi2Sr2CaCu2 08, where the surfaces are believed to be representative of the bulk, the 
densities of states at the YCBCO surfaces are shown to be qualitatively different from the bulk, and 
are sensitive to doping. The calculated spectral function agrees closely with both bulk-sensitive and 
surface-sensitive photoemission results, while the calculated density of states for optimally-doped 
YCBCO agrees closely with tunneling experiments. We find that some density of states features 
previously ascribed to competing order can be understood as band structure effects. 



I. INTRODUCTION 

The density of states (DOS) measured by tunneling 
experiments in the YBa2Cu307_5 (YBCO) family of 
high temperature superconductors is complicated. Sev- 
eral experiments on superconducting samples have mea- 
sured densities of states with multiple energy scales^"— 
Some experiments find a subgap feature^i^ while others 
do notf^'^'^'^ and all experiments near optimal doping 
find satellite features at energies larger than the gap en- 
ergy. Many of these studies find that the spectra change 
qualitatively with doping and can vary significantly at 
different points on the sample surfacei^"^ Spectral fea- 
tures have been interpreted in terms of band structure)^ 
competing order )^ and coupling to bosonic modesi^ 

In this work, we explore reasons why the YBCO 
single-particle spectrum is so complicated, particularly 
when compared to that of the related superconductor 
BisSrsCaCuaOg (BSCCO)^-^^ In BSCCO, there is a 
clear d-wave-like gap in the density of states, and it has 
been possible to reproducibly extract detailed informa- 
tion about the band structureiii^iii and superconduct- 
ing stateJ^ii^ At present, there is no consensus on how 
to interpret the tunneling DOS in YBCO. In this work, 
we focus on two specific structural differences between 
BSCCO and YBCO, namely that YBCO is a polar ma- 
terial while BSCCO is not, and that YBCO has metallic 
one-dimensional CuO chains while BSCCO does not. We 
show that the confluence of these two factors explains 
some features of the experimentally- measured DOS. In 
particular, our results suggest that some features that 
were previously thought to indicate charge-ordering ac- 
tually come from peculiarities of the YBCO band struc- 
ture. 

The polarity of the YBCO unit cell is important for 



several reasons. Unlike BSCCO, YBCO has no natural 
cleavage plane (i.e. no plane along which ionic forces van- 
ish) , making it difficult to prepare surfaces that are clean 
enough for experiments. More importantly for this work, 
there is a charge transfer between YBCO surfaces as a 
result of the electric fields generated by the polar unit 
cell. This charge transfer leads to the formation of sur- 
face states that can differ significantly from states in the 
bulk. In contrast, it is widely believed that the surface 
layers in BSCCO are representative of the bulk. 

Surface charging in polar crystals has recently be- 
come prominent in the context of LaAlOa/SrTiOs 
interfacesiiiii^ The essential idea is that, since the oppo- 
site faces of a polar unit cell have opposite charge, there 
is a potential difference between them. In a thin film, the 
potential difference between the top and bottom surfaces 
of the film is equal to the potential difference across a 
single unit cell times the number of unit cells spanning 
the film. The potential difference between the surfaces of 
the film is thus proportional to the film thickness, much 
like a parallel plate capacitor, and is typically large when 
the sample is more than a few unit cells thick. A "polar 
catastrophe" (i.e. a divergent electrostatic energy as the 
sample becomes macroscopically thick) is avoided by a 
transfer of charge between the two surfaces. This screen- 
ing charge eliminates the potential difference across the 
film, but changes the doping at the surfaces and leads to 
the formation of surface states. The existence of surface 
states in YBCO has recently been confirmed by angle 
resolved photoemission (ARPES) experimentsji^i^S 

The second aspect of YBCO that makes it distinct 
from BSCCO is the presence of layers of one-dimensional 
CuO chains, in addition to the Cu02 plane bilayers. The 
naive view is that these chains carry charge in parallel to 
the Cu02 bilayers, but have little direct impact on the 
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Cu02 layers. The chains are, therefore, generahy ignored 
in models of YBCO. For in-plane transport experiments, 
this point of view appears justified since it is possible 
to eliminate the effects of the chains (which run parallel 
to the b axis) by measuring transport in the a-axis di- 
rection; however, the chains are not easily disentangled 
from most other types of experiment. For example, c-axis 
currents (perpendicular to the planes) must pass through 
the Cu02 and CuO layers in series, so that the c-axis re- 
sistivity is dominated by the Fermi surface mismatch be- 
tween plane and chain layers.-^^ As another example, CuO 
chains have been argued to cause an anomalous vortex 
core expansion at low magnetic fields in YBCO,^^ which 
is connected to a small superconducting energy scale in 
the chains In a similar vein, we find in this work that 
the effects of the chains are subtle, but are important for 
understanding some details of the density of states in the 
Cu02 layers. 

In this work, we calculate the tunneling DOS for a 
tight-binding model of YBCO that is based on recent ex- 
perimental ARPES measurements of the band structure. 
There have, over the years, been many attempts to mea- 
sure the YBCO spectrum using ARPES but, for reasons 
discussed above, it is difficult to do reliably. Some of the 
first successful measurements were made by Schabel et 
0^24^25 .^jjp found a complicated set of bands, not all of 
which could be easily related to bands predicted by den- 
sity functional theory (DFT) calculations. Later work by 
Lu et al2^ identified an anisotropic spectrum consistent 
with the presence of CuO chains, and measurements in 
YBa2Cu408 confirmed the existence of a pair of chain 
Fermi surfaces in that material.^ Recent experiments by 
Zabolotnyy et al^ and Nakayama et all? mapped out 
the Fermi surfaces of the surface states in some detail. 
Most recently Okawa et a& have succeeded in imaging 
states in the bulk, and were able to partially map the 
Fermi surface and superconducting gap near the middle 
of the Brillouin zone. 

The goal of this work is to develop a self- 
consistent model for surface and bulk states in 
Yi_yCayBa2Cu307_5, and to use this model to under- 
stand details of the tunneling DOS. We remark that 
YBCO is sufficiently three-dimensional that it is not pos- 
sible to model surface states by considering a single iso- 
lated Cu02 layer, as is frequently done to model BSCCO. 
Instead, we develop a model for a c-axis oriented film of 
thickness Nc unit cells. 

The approach we take is phenomenological. In Sec. |IT1 
we develop a model Hamiltonian for YBCO from least- 
squares fits of a tight-binding dispersion to the surface 
states measured in Ref. [l^. We focus on optimal and 
overdoped YBCO, where we hope to avoid complications 
arising from competing phases or quantum critical points. 
The ARPES spectrum gives us the two-dimensional (2D) 
dispersion curves for the surface states; we assume that 
these dispersions are rigid, meaning that dispersions for 
bulk Cu02 layers may be obtained by shifting the sur- 
face dispersions up or down in energy. This process al- 
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FIG. 1: Model of a Ca-doped YBCO thin fihn. The film 
is Nc unit cells thick, and each unit cell contains two su- 
perconducting Cu02 and one metallic CuO layer, for a to- 
tal of = 3Nc conducting layers. The conducting layers 
are labelled i = 1,...,7V, with i — 1 corresponding to the 
chain-terminated surface, and i = N corresponding to the 
Cu02 plane-terminated surface. The Cu02 and CuO layers 
are described by two- and one-dimensional dispersions fip(k) 
and ^ic(k) respectively, where k = {kx,ky) and are coupled 
by interlayer hopping matrix elements t±p (plane-plane) and 
t±c (chain-plane). The nonconducting layers are not included 
in the band structure calculations, although the Y/Ca layers 
are included in calculations of the electrostatic potential. The 
BaO layers are electrically neutral and are ignored. Results 
in this work are shown for Nc = 10 unit cells. 



lows us to infer the structure of the bulk bands based on 
the measured surface states. We find that the inferred 
bands agree with the measured bulk Fermi surface from 
Ref. |2^ We then perform self-consistent calculations for 
the electrostatic potential in a finite-thickness film; this 
gives us the band-bending profile near the surfaces. We 
add superconductivity to the Cu02 layers by hand, us- 
ing a phenomenological relationship between the magni- 
tude of the order parameter and the charge density. In 
this way, it is found that the order parameter is smaller 
on Cu02-terminated surfaces than in the bulk. Finally, 
we calculate the superconducting DOS and normal state 
spectral functions for both surface and bulk layers. We 
discuss the results of these calculations in Sec. lIIII In this 
section, we conclude that, unlike in BSCCO, the density 
of states is sensitive to details of the band structure, and 
that these details naturally explain some of the features 
seen experimentally. 



II. METHOD 

We consider a thin film consisting of Nc unit cells, 
stacked along their c-axis, as illustrated in Fig. [TJ Each 
unit cell contains three conducting layers (two Cu02 
planes and one CuO chain layer) and three insulating 
layers (a Yi_yCaj, layer and two BaO layers). The total 
number of conducting layers is therefore N = 3Nc. Ex- 
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perimentally, it is known that YBCO cleaves at the BaO 
layer, so that the top conducting layer may be either a 
CuO chain layer or a Cu02 plane layer. To study both 
cases, we assume that the first layer (i = 1) is composed 
of CuO chains, and the last layer {i = N) is a Cu02 
plane. 

Of the three nonconducting layers, two (the BaO lay- 
ers) are nominally neutral and are not explicitly consid- 
ered in the model. The yttrium and calcium atoms are 
nominally in the and Ca^+ states and are retained 
in calculations of the long-ranged Coulomb potential. 

Of the three conducting layers, it is assumed that 
only the CuO 2 planes are intrinsically superconduct- 
ing, and that the chains are superconducting because 



of their proximity to the planes. The proximity model 
for chain superconductivity has been discussed at length 
elsewhere.—"— The essential idea is that single-electron 
hopping between the chain and plane layers leads to mix- 
ing of the electronic wavefunctions, and to induced pair- 
ing in the chain layer. 

The mean-field Hamiltonian of the system is 

H = J2clHkCk (1) 

k 

where Cl = [cj^.,,, ci_ki, • • . , cj^^^, CN-ki] and c]■^^^ is the 
creation operator for a spin-a electron in the conducting 
layer i with 2D wave vector k — {kx, ky). is given by 
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where ^ic(p) (k) is the 2D chain (plane) energy disper- 
sion in conducting layer i. The matrix elements t±c and 
t±p are for plane-chain and plane-plane hopping respec- 
tively. We remark that, in BSCCO, the plane-plane hop- 
ping matrix element has a strong k-dependence; ARPES 
experiments^^ and DFT calculations^'^ agree that the k- 
dependence in YBCO is weaker, and we have neglected 
it here as a way of reducing the number of fitting param- 
eters in our model. 

For the intralayer dispersions, we use tight-binding 
models with nearest neighbor and third-nearest neighbor 
hopping for the chains and planes respectivley: 

Cic(k) = -2icCosfc.y + (3) 

^ip(k) = —2tp[cos kx + cos ky + 2t' cos kx cos ky 

+ t"{cos 2kx + cos 2ky) 

+ 2i"'(cos 2kx cos ky + cos 2ky cos kx)] 

+ ^^ (4) 

where the values of tc, tp, t', t", t'" , and t±p are shown 
in tablelH The potential $i includes short and long range 
Coulomb interactions, as well as offset potentials for the 
plane and chain layers, and the chemical potential. The 
method of calculating $i is discussed below. 

It is known that the superconducting order parame- 
ter in the Cu02 planes has dx^-y^ + s symmetry^ as 
a result of the orthorhombic crystal structure. From 
experiments,^ the s-wave component is approximately 



15% of the d-wave component, and we assume that this 
ratio holds at both the surfaces and in the bulk. Denot- 
ing the d-wave and s-wave components in layer i by A^^ 
and Ais respectively, we have 

A - / ^{cos{kx) - cos{ky)] + Ais, i G plane . ^ 
^'•^ " \ 0, ie chain 

Note that, although the order parameter is zero in the 
chains, superconductivity is induced by the proximity ef- 
fect. We adopt the phenomenological expression for A^^ 
in the Cu02 layers 

100[ni - ninin] „ 

A^d = — meV, (6) 

^ ^min 

where rii is the charge density in plane i and nmin = 0.7. 
We expect ([6]) to be valid on the overdoped side of the 
phase diagram. At optimal doping (with planar charge 
density Up sa 0.84 in the bulk), this gives a c?-wave order 
parameter of magnitude 46 meV, which is close to values 
inferred from recent ARPES measurements.^^ Equation 
([6]) also implies that the gap vanishes when the electron 
concentration is less than nmin- The value of nmin is not 
well-known and is probably material-dependent (depend- 
ing, for example, on the level of doping-related disorder); 
the canonical form for La2-xSrxCu04 gives nmin — 0.73, 
but tunneling experiments- on Yi_j,CayBa2Cu307_5 
suggest that superconductivity is still present at this dop- 
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ing level in YBCO. Our choice seems reasonably consis- 
tent with these experiments. 

The potentials in ([3]) and Q are calculated from a 
self-consistent mean- field treatment. We write 



Ep, i E plane 
Eci * G chain 



(7) 



where is the electrostatic potential, and ec{p) includes 
the chemical potential and the energy of the chain (plane) 
tight binding orbitals. We determine €p and €c by spec- 
ifying the bulk plane and chain charge densities at opti- 
mal doping {6 = 0.08), which we take to be Up = 0.84 
and Tic — 0.48 electrons per 2D unit cell respectively. 
Note that, once Up is chosen, is set by the constraint 
of charge neutrality, given by Eq. pT|) below. Our self- 
consistent calculations then find ec — = 3.176 eV at 
optimal doping. We assume that the tight binding or- 
bitals are not modified by doping, so that ec — fp is held 
constant througout this work. In our calculations chemi- 
cal doping modifies the band structure only through the 
electrostatic potential (f>i. 

The Coulomb potential (j>i is then calculated self- 
consistently within the Hartree approximation, under 
the assumption that the charge is uniformly distributed 
within each layer. For the planar geometry shown in 

Fig.m 



U^, i £ plane 
i e chain 



(8) 



where the Yi_yCaj, layers are implicitly included in the 
sum over j and 



(9) 



In ([8]), U is the intraorbital Coulomb potential for the 
Cu02 layers, which we take to be 4 eV. The total 2D 
charge density in layer i is ai = — rii, where is the 
charge density of the ionic cores. Here, rii is measured 
relative to the Cu'^"'" and O^^ states, so that 




j G CuOi_5 chain 
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(10) 



All charge densities are in units of e/a^, where ao ~ 4 A 
is the 2D lattice constant. Electrical neutrality requires 
that 



2np^2 + 2S ~y. 



(11) 



In Eq. ([5]), Zi is the 0-coordinate of layer i, in units of 
the c-axis lattice constant dz ~ 12 A. Within a unit cell, 
the layers are at z = (chain), z — 0.354 (first plane), 
z — 0.5 (Y layer), and z — 0.646 (second plane);^ The 
weak doping dependence of these values is ignored here. 
The dielectric constant e in ([SJ is not well known, but is 



TABLE I: Parameters for the tiglit binding model for YBCO. 



Parameter 


Value 




105 meV 


t' 


-0.277 


t" 


0.234 


t'" 


-0.042 


tc 


500 meV 


t±p 


61 meV 




l.lt±p 



believed to be around e = 20, which is the value taken 
here. 

For a given potential $j, the charge density in layer 
i is found from the eigenvalues -E^.k and eigenstates 
^Q,k(i,cr) of Hk via 



2N 



■EEl*".k(*,t)lV(i^a,k) 



(12) 



k Q = l 



where f{x) is the Fermi function, the factor of 2 is for 
spin, and Nk is the number of k-points in the sum. The 
band index a ranges from 1 to 2^^^ because the number 
of bands in the film is equal to twice (including spin) 
the number of conducting layers. The updated charge 
density is used to re-calculate which is then used in 
the next iteration for n^. The iterations proceed until the 
difference between rii in two consecutive iterations is less 
than 10~^. To reduce the computational workload, is 
calculated in the normal state. 

We finish this section with a brief discussion of the fit- 
ting procedure used to get the model parameters shown 
in Table HI As discussed above, we fit the 2D disper- 
sions for the Cu02 layers to the surface states mea- 
sured in Ref. [l9l In order to avoid complications from 
the chains, we fit the energy spectrum for an isolated 
Cu02 bilayer to the measured bands in regions of the 
Brillouin zone far from the chain Fermi surface. The 
model bilayer has bonding and antibonding bands, with 
energies ^^(k) = ^p(k) ± t±p, which allows us to de- 
termine tp, t\ t" , t'" , and t_\_p. One potential diffi- 
culty with this fitting process is that it assumes that the 
electrostatic potential is the same in the top two Cu02 
layers, meaning that we attribute the experimental bi- 
layer splitting entirely to t^p. If we allow for a poten- 
tial difference A(/) between the Cu02 layers making up 
the bilayer (due to band bending at the surfaces), then 

^±(k) = ^p(k) ± ^(A(/i/2)2 + Self-consistent cal- 
culations reported in the next section suggest A(/) 10 
meV while the measured t_Lp « 60 meV; it follows that 
A0 modifies the band energies by less than 1% and can 
safely be neglected. 

There are, at present, no reliable measurements of the 
chain band structure. We therefore assume that strong 
correlations are not significant in the CuO chains (which 
are roughly quarter- filled) , and that the chain bandwidth 



5 




k k 

X X 

FIG. 2: (Color online) Comparison of the self-consistently 
calculated spectral function at to the experimentally mea- 
sured Fermi surface (circles). Results are for (a) surface states 
(data from Ref. [l9h and (b) bulk states (data from Ref . l28h . 

of 2 eV from DFT is reasonable)^ This gives the intra- 
chain hopping matrix element tc = 500 meV. The hardest 
parameter to establish is ij_c, the chain-plane coupling. 
We show results for multiple values of t±c in Sec lIIIB| 
and find that the experimental density of states is rea- 
sonably well fit for t±c = for YBa2Cu3 06.92- 

We show the results of our fitting procedure in Fig. [2] 
In this figure, we compare calculated surface and bulk 
spectral functions at the Fermi energy with the surface 
and bulk spectral functions measured by ARPES. The 
calculations are in good agreement with the experiments, 
which suggests that the assumptions made in developing 
our model are reasonable. 



III. RESULTS AND DISCUSSION 

A. Self-consistent potential, charge density, and 
superconducting order parameter 

Figure [3] shows the results of self-consistent calcula- 
tions for the charge density and electrostatic potential. 
We take the specific case of Yi_j,CayBa2Cu306.92, which 
corresponds to optimal doping when y = and to over- 
doping when y > 0. The self-consistently determined 
charge density is shown in Fig. ^a) as a function of 
layer index i. (Recall that the layer indices i = 1 and 
i = 30 label the CuO chain-terminated and the Cu02 
plane-terminated surfaces respectively.) As discussed in 
the introduction. Fig. ^a) shows that there is charge 
transfer from the plane-terminated surface to the chain- 
terminated surface. This charge screens the electric field 
produced by the polar unit cells, so that the electric po- 
tential is constant across the thin film except near the 
surfaces. This is shown in Fig. [3jc) and (d). In these 
figures, we have plotted the difference 

A$,p(c) (13) 

between the potential in layer i and the potential for a 
plane (chain) in the bulk in order to make the comparison 




layer index i layer index i 

FIG. 3: (Color online) Self-consistent solutions for the charge 
density and electrostatic potential in Yi_yCayBa2Cu3 06.92. 
Charge density is shown (a) for all CuO chain and Cu02 plane 
layers in the thin film as a function of the layer index i, and 
(b) at the chain {i — 1) and plane (i = 30) surfaces as a 
function of y. The potential difference A$ip(c) = — *&p(c), 
between the potential in layer i and the potential in the bulk 
planes (chains), is shown in (c) and (d). 



between different y values simpler. A$jp(c) is nonzero 
within a screening length of the surfaces and is posi- 
tive (negative) at the plane (chain) surface. The charge 
density on the plane (chain) surface is correspondingly 
smaller (larger) than in the bulk as shown in Fig.[3Ja). 

We note that A$ic, the potential shift at the chain- 
terminated surface, is roughly seven times larger than 
A^30pj the potential shift at the plane-terminated sur- 
face. This ratio is approximately the same as the ratio 
of the plane and chain densities of states (in the nor- 
mal state) ; because of the small density of states in the 
chains, a large chemical potential shift is required at 
the chain-terminated surface to accommodate the charge 
transferred from the plane-terminated surface. 

Figure [3] also shows the effect of Ca substitution. De- 
spite the proximity of the Ca ions to the CuO 2 layers, rii 
changes by roughly the same amount in both the Cu02 
and CuO layers. The one notable exception to this is 
at the chain surface, where the electron concentration 
changes by roughly half as much as in the bulk [Fig. ^h)] . 

Figure 2] shows the d-wave superconducting order pa- 
rameter Aid as a function of layer index i in the Cu02 
plane layers, and for different levels of Ca doping. Be- 
cause Aid is calculated phenomenologically from Eq. ^ , 
it follows Tii. Thus, Aid is larger than in the bulk near 
the chain surface and smaller than in the bulk near the 
plane surface. In our calculations. Aid actually vanishes 
at the Cu02 surface layer (z = 30), although there is 
a spectral gap due to proximity coupling to the subsur- 
face layers. Note that recent ARPES experiments,—!^ 
found that the surface is nonsuperconducting, but tun- 
neling experiments'':^^ found a clear superconducting gap 
at the surface. We will show in the next section that the 
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FIG. 4: (Color online) Superconducting order parameter as a 
function of layer index for Yi_yCayBa2Cu3 06.92- The figure 
shows the d-wave component of the order parameter, which 
is nonzero in the Cu02 layers only. The magnitude of A^^i 
is given by Eq. ((6|. Inset: Superconducting order parameter 
near the plane surface (i — 29) as a function of Ca doping. 
Note that Asod = 0. 




(0 (meV) 



FIG. 5: (Color online) Surface and bulk density of states for 
Yi_yCayBa2Cu3 06.92. Results are shown for y = 0.0 [(a) and 
(b)] and for y = 0.05 [(c) and (d)]. DOS is shown for Cu02 
planes [(a) and (c)] and [(b) and (d)] CuO chains. Model 
parameters are given in Tabled 



size of the superconducting gap at the Cu02 surface is 
sensitive to how the surface is prepared. 

B. Density of states and spectral function 

The main results reported in this work are for the layer- 
dependent density of states and spectral function. The 
density of states at energy uj in layer i is given by 




where TVk is the total number of k-points, and Ai{'k,uj) 
is spectral function in layer i, 

2N 

+ |*a,k(«,i)P'5(c^ + i;„,k)]. (15) 

Note that the eigenstates ^Q,,k(*, cr) and eigenenergies 
Ea,ii describe electrons {a =t) or holes {a =1) depending 
on the spin index. In this work, all results for Ai(k,uj) 
are shown for the normal state at uj = ep, while results 
for pi{uj) are shown in the superconducting state. 

Results for Yi_j^CayBa2Cu3 06.92 with y = and 
y = 0.05 are shown in Fig. [S] This figure shows that 
the density of states at the surfaces is very different from 
in the bulk, for both the Cu02 planes and CuO chains. 
Notably, the d-wave gap is significantly smaller at the 
surfaces than in the bulk. In the Cu02 layers, this re- 
flects the suppression of A^c; near the surface due to the 
decreased electron density in the surface layers. This is 



consistent with bulk-sensitive ARPES experiments that 
found a gap of order 42 meV,'^® and surface-sensitive tun- 
neling experiments, which consistently find that the gap 
is ^ 25 meVii"— In the chain layers, the reduced gap 
at the surface reflects the reduced proximity coupling at 
the chain surface, in part because the chain has only one 
nearest neighbor Cu02 plane. Our calculation shows that 
the DOS obtains its bulk value within a few layers of ei- 
ther surface. 

We note that, relative to the conventional model of a 
single-layer d-wave superconductor, the DOS in Fig. [5] 
shows a lot of structure. In particular, the Cu02 surface 
(for y = 0) has a main gap of about 10 meV, whose co- 
herence peaks are split into pairs of closely-spaced peaks, 
and a satellite "shoulder" at about 20 meV. Shoulder fea- 
tures have been commonly observed in tunneling experi- 
ments, and have been attributed to pairing,*^ and to com- 
peting order. ^ In our calculations, this structure comes 
from the interplay of pairing and band structure effects, 
namely the mixing of chain and plane states resulting in 
an orthorhombic distortion of the Fermi surfaces. This 
is illustrated in Fig. [51 which shows the layer-resolved 
spectral function Aiik^ep) at the Fermi energy. We see 
from this figure that chain-plane coupling strongly dis- 
torts the Fermi surfaces in the (tt, 0) region of the Bril- 
louin zone. This distortion is particularly important for 
the Cu02 layers, because there is a van Hove singularity 
near the (tt, 0) point in the undistorted spectrum. The 
DOS is therefore sensitive to small changes in the Fermi 
surface shape, caused either by chain-plane coupling, by 
doping, or by band-bending (changes in the electrostatic 
potential) near the surfaces. Thus, the addition of 5% 
Ca changes the electron density by only ~ 0.02 electrons 
per Cu02 plaquette but qualitatively changes the shape 
of the coherence peaks [compare Fig. [5ja) and (c)]; in 
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FIG. 6: Spectral function at the surfaces and in the bulk of 
Yi-j,CayBa2Cu3 06.92. Columns show A(k, ef) for y = for 
(a) planes and (b) chains, and for y = 0.05 for (c) planes and 
(d) chains. Rows (top to bottom) show layers i = 30, 29, 27, 
15 (planes) and i — 1, 4, 7, 13 (chains); thus, surface states 
are shown in the top row, states deep in the bulk are shown 
in the bottom row. 



Fig. ini we see that the main effect on the Cu02 surface 
layer {i = 30) is indeed near the (tt, 0) point, where the 
spectral weight is reduced by Ca doping. 

In Fig. [HI the intensity of A,j(k, w) strongly depends 
on the amount of hybridization between chain and plane 
states. The Fermi surface at the CuO-terminated surface 
{i = 1) is relatively undistorted, indicating weak chain- 
plane coupling. By contrast, the chain Fermi surface in 
the bulk is more strongly hybridized to the Cu02 plane 
states and is correspondingly washed out. This is a pos- 
sible reason that the chain Fermi surface was imaged in 
surface-sensitive ARPES measurements^—"^ but not in 
bulk-sensitive measurements 

We remark that the amount of hybridization between 
chains and planes is strongly k-dependent, even though 
the matrix element t±c is independent of k. This is be- 
cause the hybridization at each k depends on the energy 
difference between the chain and plane bands at that 
value of k. Because these bands have different symme- 
tries, the energy difference (and thus the hybridization) 
is a strong function of k. It follows that the amount of 
hybridization at ej? is a strong function of the relative 
positions of chain and plane Fermi surfaces, with the hy- 
bridization being largest where the Fermi surfaces cross. 
This explains the difference between the surface {i — 1) 
and bulk (j — 13) spectral functions shown in Fig. [S) A 
more extensive discussion of this point can be found in 
Ref. M 

The degree to which the chain and plane Fermi sur- 
faces hybridize determines the size of the induced gap 
in the chain layers. The small gap at the CuO surface 
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FIG. 7: Effects of an adsorbed surface layer on the surface 
states of YBa2Cu3 06.92. Columns are for (a) plane and (b) 
chain DOS, (c) plane and (d) chain spectral function. The 
adsorbate layer has a 2D charge density of +ena for the Cu02 
surface and —eus for the CuO surface. Rows correspond to 
Us = 0.0, 0.1, 0.15, and 0.2 (top to bottom). 



(Fig. [5]) is thus due to the weak hybridization of the sur- 
faces chains with the underlying Cu02 planes. In the 
bulk, the chain DOS shown in Fig. [S] has a main gap 
of about 20 meV, and a small gap of about 4 meV. The 
small gap originates from sections of the chain Fermi sur- 
face that are only weakly coupled to the planes (namely, 
^a; ^ 1) while the large gap comes from sections of the 
Fermi surface that are strongly coupled to the planes 
{kx ^ 1). The small gap was discussed previously as 
a possible source for the subgap structure measured in 
some tunneling experiments;^ 

Having discussed general features of the DOS and spec- 
tral function, we now discuss how these are affected by 
changes in specific model parameters. First, we allow for 
the possibility that the surface layers of the YBCO thin 
film are doped by the adsorption of atoms or molecules 
onto the surfaces. Adsorption happens naturally, for ex- 
ample, when YBCO is exposed to air, and deliberate 
potassium adsorption has been used to control the elec- 
tron concentration in Cu02 surface states. '^^ In this work, 
we are interested in the possibility that the adsorbate 
layers partially screen the electric fields at the YBCO 
surfaces. We model this by assuming that the adsorbed 
layer has an average charge density of +ens {—ens) at 
the Cu02 (CuO) surface, where Ug > 0. In our self- 
consistent calculations, the effect of the adsorbate charge 
is to reduce the charge transfer between the chain and 
plane surfaces. 

Figure [7] shows the effect of on the density of states 
and spectral function. We see that even a relatively small 
adsorbate charge density has a significant effect on both 
Pi{uj) and Ai(k,uj) at the surfaces. In particular, pi{u}) 
and Ai{\i,uj) at the surfaces are increasingly similar to 
the bulk as rig increases. This follows directly from the 
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FIG. 8: Effects of chain-plane coupling on the surface states 
of YB2Cu3 06.92- Columns are for (a) plane and (b) chain 
DOS and (c) plane and (d) chain spectral function. The ad- 
sorbate charge density is Us = 0.15. Rows correspond to 
t±c = 0.0, 0.8t±p, l.lt±p, 1.3t±p (top to bottom). 



fact that the charge densities at the surfaces approach 
their bulk values as increases. In the Cu02 layers, 
this results in a larger Aid from Eq. ([6]) , while in the CuO 
layers, this results in an increased hybridization between 
the surface chains and the adjacent Cu02 plane. 

The Cu02 spectrum for Ug = 0.15 shown in Fig. [7| 
is consistent with existing tunneling experiments on 
optimally-doped YBCO. For this case, we obtain co- 
herence peaks at the Cu02 surface at ±20 meV and a 
satellite peak at « —30 meV, in approximate agreement 
with Refs. 01111 (when comparing with experiments, re- 
call that a peak at negative energy in the DOS corre- 
sponds to a peak at positive voltage bias in a tunneling 
experiment). We note that having a charged adsorbate 
layer is not the only way to obtain agreement with exper- 
iments; a different value for the dielectric constant, for 
example, will affect the surface charge density, and conse- 
quently the DOS. Whether or not our model is correct in 
all details, it demonstrates that existing tunneling mea- 
surements on optimally-doped YBCO can be explained 
within a band picture. 

Throughout this work, we have assumed that the 
plane-chain coupling parameter is t±c — l-l^±p = 
67 meV. This is the hardest of the model parameters 
to establish, and was chosen because it gives reasonable 
results for pi{uj) and ^^(k, cj). Figure [8] shows the ef- 
fect of varying t±c on the density of states and spectral 
function. Not surprisingly, chain-plane coupling has lit- 
tle effect on the DOS at the Cu02 surface, apart from a 
shift of a weak negative-energy satellite peak away from 
the Fermi energy with increasing t±c- In contrast, the 
chain surface is strongly influenced by coupling to the 
Cu02 layer; the CuO chains are metallic when t±c — 0, 
and an induced gap appears when t±c is nonzero. The 
induced gap grows approximately linearly with t±c and 




CO (meV) 

FIG. 9: Effects of Ca doping on the density of states of 
Yi_yCaj,Ba2Cu306.92 at the surfaces. Columns are for (a) 
plane and (b) chain surfaces with Us = 0.0, and (c) plane 
and (d) chain surfaces with Us — 0.15. Rows correspond to 
y = 0.0, 0.05, 0.10, 0.2 (top to bottom). Arrows in (d) indi- 
cate the locations of satellite features discussed in the text. 




FIG. 10: Effect of Ca doping on the spectral function at the 
surfaces of Yi_yCaj;Ba2Cu3 06.92. Columns are for Us — 0.0 
for (a) planes and (b) chains, and for Us — 0.15 for (c) planes 
and (d) chains. Rows are for y = 0.0, 0.05, 0.10, 0.2 (top to 
bottom). 



has both a subgap and a main gap, as discussed earlier. 
The spectral function also changes with increasing t±c, 
becoming increasingly distorted near (tt, 0). 

Finally, we look more closely at the effects of Ca 
doping. In Ngai et al^ samples with up to 20% Ca 
substitution for Y were studied, while in Yeh et al^, 
30% Ca-doped samples were studied. In Fig. ^ we 
show a series of calculations for the surface DOS of a 
Yi_yCaj,Ba2Cu306.92 thin film. Results are shown for 
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both Us = and rig = 0.15. The corresponding spectral 
functions are shown in Fig. [TUl 

As discussed above, the DOS exhibits subgap, main 
gap, and satellite features. In Ngai et al, the three fea- 
tures were attributed to superconductivity on different 
regions of the Fermi surface. This is generally consistent 
with our findings here although we have found it diffi- 
cult to attribute some satellite peaks to specific Fermi 
surface elements. Experimentally, the ratio between the 
satellite, main gap, and subgap energies was found to be 
approximately constant as a function of Ca doping, and 
this was argued to show that there is a common pairing 
mechanism for all three features. In our model the en- 
ergies of the satellite and subgap features in Fig. [9] also 
scale with the main gap. 

In some cases, our calculations reproduce the detailed 
structure of the experimental spectra. For y — 0.05, Ngai 
et al showed two types of spectrum. The first is remark- 
ably similar to that shown for y — 0.05 in Fig. IHJc), 
having pronounced coherence peaks and satellite features 
resembling shoulders. Spectra of this type are also mea- 
sured in optimally-doped {y = 0) samplesi^i^i^i^ The good 
agreement between our calculations and the measured 
spectra suggests that our model captures the essential 
physics of the surface states in optimally-doped YBCO. 

The second type of measured spectrum qualitatively 
resembles that shown for y — 0.0 in Fig. [SJd), having 
weak coherence peaks, a subgap feature, and satellite 
peaks (indicated by arrows in Fig. [S]). However, there 
is a discrepancy between our calculations and the exper- 
iments; at higher Ca doping levels, the measured spectra 
continue to exhibit three sets of peaks while the satellite 
peaks in our calculations become less prominent. Given 
the sensitivity of the DOS to small changes in the model 
parameters, it is plausible that this discrepancy can be 
corrected by small changes to the model. It is also possi- 
ble that extrinsic effects not considered here, for example 
tunneling matrix elements that emphasize the (tt, 0) and 
(0, tt) regions of the Brillouin zone^i^ could increase the 
prominence of the satellite features in tunneling experi- 
ments. However, we also cannot rule out the possibility 
that our simple model for Ca substitution is overly naive. 

Another feature of the DOS that is not explained by 
our model is the large residual DOS measured experi- 
mentally; the DOS is never seen to vanish at ep in the 
superconducting state, and is often 50% of the normal 
state DOS. It is not clear whether the residual DOS is 
intrinsic (for example, due to pair breaking at the sur- 
face) or extrinsic (coming from surface states in an ad- 
sorbate layer). It is possible that a full description of the 
YBCO surface states will require a proper accounting of 
this residual DOS. 



IV. CONCLUSIONS 

We have studied the surface and bulk states of 
Yi_yCayBa2Cu3 06.92 with < y < 0.2 within a tight- 



binding model. The model parameters for the Cu02 
planes are extracted from photoemission experiments, 
and self-consistent calculations are used to relate the sur- 
face and bulk states. We have calculated the density of 
states Pi{uj) and spectral function Ai(k,eF) as functions 
of adsorbed surface charge density, chain-plane coupling, 
and Ca doping. Our main findings are that 

• our model produces results which are in simultane- 
ous agreement with surface and bulk ARPES mea- 
surements. This supports two key assumptions of 
the model, that the surface and bulk bands are con- 
nected by a simple chemical potential shift and that 
the Coulomb potential can be treated in a planar 
approximation. 



• the DOS in optimally-doped YBCO can be quan- 
titatively explained by our model. In particular, 
shoulders measured in the density of states that 
were previously attributed to pairing or to compet- 
ing order are found to be band structure effects. 



• the superconducting DOS is sensitive to small 
changes in the model parameters. This suggests 
that, in materials with CuO chains, small changes 
in doping of the surface states can have a qualita- 
tive effect on the density of states, purely as a result 
of changes to the band structure. This should be 
contrasted with BSCCO where it is believed that 
the doping dependence of the DOS is primarily due 
to strong correlations. 



• we can understand some features of Ca-doped 
YBCO; for example, our calculations find densities 
of states at the chain surface with subgap, main 
gap, and satellite features similar to experiments. 
However, we have not understood the dependence 
of this spectrum on Ca concentration. 
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